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INTRODUCTION 


The  principal  obstacle  to  the  construction  of  accurate  global  tide  models 
Is  the  lack  of  observations  In  mld-occan.  There  Is  currently  a  wealth  of  tide 
gauge  data  from  the  shores  and  estuaries  of  most  developed  countries  and  a  few 
scattered  Islands  In  ald-occan.  These  data  have  occasionally  been  supplemented 
by  data  from  bottom  pressure  gauges  emplaced  for  relatively  short  periods  of 
time  at  a  very  few  locations.  Recently,  land-based  tidal  gravity  data  have 
also  been  used  In  tidal  modeling  (Kuo  and  Jachens  [1977]).  But  to  model  deep 
ocean  tides  from  these  essentially  land-based  data  Is  analogous  to  modeling 
a  hurricane  In  mid-ocean  by  observing  swell  at  the  shore.  Only  the  grossest 
features  can  be  depicted  and  many  assumptions  must  be  made  which  limit  the 
accuracy  of  such  models.  Tills  Is  evident  from  a  review  of  present  tide  model¬ 
ing  techniques. 

Many  Investigators  (Pekerls  and  Accad  [1969],  Zahol  [1970],  Hcndershott 
[1966],  Estes  [1977])  have  developed  global  tide  models  by  numerical 
integration  of  LaPlace's  tidal  eruatlons,  but  these  models  require  additional 
models  (or  assumptions)  of  related  physical  parameters  such  as  coastline  geo¬ 
metry,  bottom  topography,  energy  dissipation,  earth  tides,  scl f-gravitat ion , 
and  crustal  loading.  Their  sensitivity  to  uncertainty  in  these  parameters  can 
be  dramatic.  Pekerls  and  Accad  [1969]  find  that  small  changes  in  dissipation 
or  coastline  parameters  can  cause  amplitudes  to  change  by  a  factor  of  3  at J 
aaphldromca  to  rotate  In  the  wrong  direction.  As  a  result,  Kuo  and  Jachens 
[1977]  among  others,  find  that  such  models  generally  give  poor  agreement  with 
t  ld.il  observations  on  islands  in  mid-ocean. 

In  an  alternative  approach,  Kuc  and  Jachens  [1977]  recently  used  tidal 
gravity  ohservat  ions  to  model  open  ocean  tides  in  limited  regions  of  the 
world's  oceans.  This  approach  models  only  the  relatively  long  wavelength 
features  of  the  ocean  tides  (typically  1800  km  or  more  In  scale  length) 
depending  on  the  number  and  distribution  of  land  based  tidal  gravity  stations. 
Since  it  is  an  inversion  technique,  it  docs  not  yield  unique  solutions  and 
must  be  constrained  by  bottom  pressure  observations,  and  it  is  sensitive  to 
the  locat  ion  of  the  tidal  gravity  stations.  An  important  assumption  of  this 
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technique  is  that  the  tide  outside  the  region  to  be  modeled  is  accurately 
known,  and  thus  this  technique  is  dependent  on  other  tidal  models  for  the 
rest  of  the  globe.  Other  regional  tide  models  using  bottom  pressure  gauges 
(Hunk,  et  al  [l97o]),  require  unrealistic  smoothing  of  the  bottom  topography 
(Hendershott ,  [l97l]). 

Many  of  the  above  mentioned  problems  may  be  solved  by  the  use  of  direct 
satellite  altimetry  measurements  (Zetler  and  Maul  J1971)).  Since  the  satellite 
altitude  measurement  is  referenced  to  the  earth's  center  of  mass,  such  measure¬ 
ments  yield  the  true  elevation  of  sea  surface  relative  to  this  point,  and  are 
free  of  assumptions  about  bottom  topography,  earth  tides,  crustal  loading, 
coastlines,  and  the  distribution  of  terrestrial  measurement  stations.  In 
addition,  satellite  measurements  can  be  made  at  high  spatial  density  over  all 
the  world's  oceans  in  a  short  time.  Figure  1  shows  the  distribution  of  CEOS-3 
satellite  altimeter  measurement s  resulting  from  the  first  year  of  operation. 

This  satellite  was  Just  recently  turned  off  after  more  than  A  years  of  operation. 
The  Scasat  satellite,  which  was  launched  in  June  of  1978,  produced  altimeter 
data  which  supplements  the  CEOS-3  data  distribution.  Even  though  this  satel¬ 
lite  operated  only  99  days,  altimeter  data  to  a  spatial  density  of  50  km,  from 
72°N  to  72°S  latitude  was  produced  uniformly  over  the  globe.  This  data  serves 
to  fill  in  the  gaps  in  the  CEOS-3  data  distribution. 

While  the  altimetric  technique  showed  great  potential,  it  also  had  formid¬ 
able  practical  problems.  The  simulation  study  by  Zetler  and  Maul  1 1971], 
showed  that  amplitudes  and  pharos  of  tidal  components  could  be  recovered  in 
the  presence  of  realistic  satellite  orbit  errors,  but  Maul  and  Yanaway  [1978] 
found  that  actual  orbit  errors  were  too  large  and  not  random  enough  to,  permit 
successful  tide  parameter  recovery.  They  concluded  that  1  m  orbit  accuracy 
must  be  achieved  before  tidal  recovery  from  altimetry  is  possible.  Won  and 
Miller  [1978]  attempted  to  mitigate  the  effects  of  orbit  error  by  simultaneously 
.Ivin,;  for  gcoid,  tides,  and  bias  type  orbit  errors  from  altimetry  in  the 
hortb  Atlantic.  Although  they  achieved  good  agreement  with  independent  geoid 
models  and  reduced  tidal  amplitude  residuals  to  less  than  a  meter,  the  agree¬ 
ment  with  known  tide  parameters,  especially  the  M^  tide  in  the  North  Atlantic 
area,  was  disappointing.  They  cited  insufficient  altimeter  data,  excessive 
orbit  errors,  and  excessive  altimeter  data  noise  as  the  principal  reasons  for 
the  poor  tidal  parameter  recovery. 
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Recently,  however.  Brown  and  Lo  1 1978)  were  successful  In  the  determination 
of  tidal  parameters  In  the  Culf  of  Alaska,  from  analysis  of  CEOS-3  altimeter 
data,  by  using  a  unique  method  for  removing  orbit  errors.  Each  pass  of  alti¬ 
meter  data  was  first  converted  to  sea  surface  heights,  by  using  an  initial  (albeit 
erroneous)  estimate  of  satellite  orbit  ephemeris.  These  sea  surface  height 
data  are  then  fitted,  pass  by  pass,  to  a  model  geold  height  surface.  The  fit 
Is  effected  by  varying  the  satellite  orbit  parameters.  Since  each  pass  is 
fitted  Independently  to  the  same  surface,  all  long  wavelength  orbit  errors 
are  removed  while  preserving  short  wavelength  tide  signals.  After  orbit 
errors  are  removed,  an  approach  similar  to  that  of  Maul  and  Yanavay  (1978) 

Is  used  to  define  tidal  phases  and  amplitudes.  Tidal  phases  are  defined  by 
analysis  of  sea  height  differences  at  intersections  (crossovers)  of  satellite 
subtracks,  and  a  least  squares  harmonic  analysis  is  conducted  to  solve  for 
amplitudes  and  phases  of  the  tidal  components.  Estimates  of  the  tidal  phase 
and  amplitude  from  6  locations  in  this  region  were  generated  from  a  preliminary 
collection  of  data.  Results  at  these  locations,  representative  of  different 
bathymetric  regimes,  agree  in  general  with  the  bottom  pressure  gauge  measure¬ 
ments  of  Rapatr,  et  al  (1977)  in  the  Culf  of  Alaska  seamount  province,  and 
confirm  the  existence  of  4  rcter  tidal  range  in  the  Culf  of  Alaska.  In  a 
separate  test,  recovery  of  M.,  tide  parameters  at  the  location  of  the  MODE  3 
bottom  pressure  gauge  (Zetler,  et  al  (1975))  showed  agreement  to  within  26  cm 
in  amplitude  and  30  degrees  In  phase. 

?h 1  .  approach  is  extended  in  this  study  to  26  local  Ions  in  the  Northeast 
i'aclfiC  i'ci'.m.  While  it  is  sensitive  to  the  altimeter  data  accuracy  and  the  distri¬ 
bution  of  crossovers  In  space  and  time,  it  Is  much  less  sensitive  to  satellite 
orbit  errors  than  the  techniques  of  Won  and  Miller,  or  Maul  and  Yanaway.  In 
addition.  In  its  reliance  on  crossover  comparisons.  It  Is  immune  to  the  effects 
of  sharp  gradients  In  the  geold  or  bottom  topography.  Since  it  generates 
strictly  local  quasl-polnt  solutions.  It  can  be  done  quickly  and  efficiently, 
without  the  need  for  modeling  effects  In  a  larger  region.  It  can  also  create 
..lohal  tide  naps  to  a  very  fine  scale  by  doing  separate  local  point  solutions 
it  ;h<  required  spatial  separation,  and  is  not  inherently  limited  to  regional 
r  lotal  application.  The  objectives  of  this  study  are  to: 


•  develop  a  technique  for  low-cost  survey  of  ocean  tide  phenomena 
of  long  and  short  wavelength  In  mid-ocean,  completely  independent 
from  models  of  bathymetry,  coast  lines,  and  earth  lines 

•  produce  a  detailed  map  of  the  co-tidal  and  co-range  lines  of  the 
major  tidal  components  (M,,  Kj,  S,, ,  0^,  and  N,)  In  the  Northeast 
Pacific  Ocean 

•  provide  a  comparative  analysis  of  altlmetrlc  tide  models  and 
existing  side  models  in  the  Northeast  Pacific  for  validation  of 
these  models. 

In  the  sections  which  follow,  the  state  of  the  art  in  tide  modeling  in 
the  Northeast  Pacific  is  reviewed,  followed  by  a  discussion  of  data  processing 
performed  to  date  on  the  CEOS-3  altimeter  data.  The  locations  of  potential 
tide  solutions  permitted  by  the  distribution  of  CEOS-3  and  SEASAT  satellite 
altimeter  data  are  identified,  and  preliminary  solutions  of  tide  parameters 
are  discussed.  Finally,  plans  for  completing  the  validation  of  tide  models 
In  the  Northeast  Pacific  using  satellite  altimetry  arc  presented. 


TIDE  DETERMINATION  IN  THE  NORTHEAST  PACIFIC 


To  provide  a  reference  for  Judging  altlmetrlc  solutions  for  tides,  we 
review  here  the  results  of  previous  attempts  to  model  tides  in  the  North¬ 
east  Pacific  Ocean.  Rather  than  compare  the  models  everywhere,  a  single 
comon  geographic  location  is  selected.  This  point  is  chosen  in  deep  ocean, 
at  about  215°W  longitude  and  34°N  latitude,  because  bottom  pressure  guage 
measurements  were  available  for  this  location  (Irish,  et  al  (1971]).  This 
Is  the  JOSIE  II  station.  A  single  point  comparison  also  better  serves  to 
Illustrate  the  uncertainty  In  deep  ocean  tide  models,  rather  than  their 
simllarlt ies. 

Most  investigators  model  the  component  of  the  tide.  This  is  the 
largest  component  of  the  tides,  but  unfortunately,  the  least  well  determined 
Figure  2  illustrates  this  uncertainty.  amplitudes  (indicated  by  the  arrow 

length)  and  phase  angles  (Indicated  by  the  arrow  direction)  range  from  13  to 
57  cm  and  are  scattered  over  a  300°  arc  of  phase  angles.  This  situation  is 
due  to  two  factors,  the  near  resonance  of  the  semidiurnal  tide  (M^,  ,  N?.  etc.) 

in  this  part  of  the  ocean  ( Hendershot t ,  [1973]),  and  the  indirect  tide  determination 
techniques  used  by  most  investigators.  Since  most  investigators  use  tide 
measurements  which  are  rather  far  removed  from  the  deep  ocean,  many  assumptions 
and  auxiliary  physical  models  are  necessary  in  the  calculation  of  the  tide  in 
deep  water.  Small  changes  in  these  models,  such  as  coast  line  configuration, 
ocean  depth,  dissipation,  earth  tides,  ocean  bottom  loading,  and  sel f -at t ract ion , 
can  lead  to  tidal  amplitudes  that  fluctuate  by  3  meters  or  more,  and  amphldrones 
that  shift  location  or  even  rotate  in  the  wrong  direction. 

It  is  instructive  to  review  the  techniques  used  In  these  tide  determina- 
t‘  '-i  .  Proceeding  in  more  or  less  chronological  order,  the  earliest  techniques 
ar  t  <•  sent -enp  lr  le  a  1  techniques  of  Harris  11904]  and  Dietrich  (1944],  Only 
ph.i •.«  angles  (co-tidal  lines),  not  amplitudes,  were  determined.  Observed  tide 
phases  at  the  coasts  and  at  some  islands  were  used  to  determine  the  co-tidal 
lines  of  awphldroBiic  systems.  These  lines  were  extrapolated  to  determine 
the  locations  of  the  anphidrome. 
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Bogdanov  et  ul  (1961a, bl,  (1964)  used  the  harmonic  constants  of  the  tide 
at  coastal  und  Island  stations  to  determine  the  Instantaneous  relief  of  the 
ocean  surface  (lsohypse  method).  Co-tidal  lines  were  determined  by  comparing 
lsohypse  charts  at  appropriate  times.  Where  data  were  lacking,  as  on  the 
southern  boundary  of  the  world's  oceans.  Interpolated  values  were  used. 

Hendershott  (1966)  used  the  principles  of  hydrodynamics  (numerical  solu¬ 
tion  of  Laplace's  tidal  equation)  in  a  global  solution  for  tides  in  Idealized 
ocean  basins.  Observed  tides  at  coastal  stations  were  used  as  boundary  eondi- 
t  ions. 

Tlron  et  al  [ 1967 1  also  used  a  global  numerical  solution  a  la  Kendershott, 
but  included  island  data  as  well.  Dissipation  was  parameterized  by  specifying 
coastal  elevations,  anu  earth  tides  were  ignored. 

The  tide  determinations  by  Munk  et  al  (1970)  were  local  in  nature,  confined 
to  the  area  along  the  California  coast  between  San  Francisco  and  San  Diego,  and 
extending  westward  halfway  to  Hawaii.  Rather  than  relying  on  coastal  measure¬ 
ments,  which  they  regarded  as  anomalistic  for  deep  ocean  modeling,  they  used 
bottom  pressure  gauge  measurements  at  eight  stations  in  deep  water  beyond  the 
continental  shelf,  the  model  is  a  normal  mode  representation  of  the  coastal 
tide,  without  dissipation,  but  accounting  for  the  earth  tide  due  to  astronomical 
forces.  Their  model  predicted  an  amphidrome  at  about  27°?!,  135°W.  To  test 
this  prediction,  Irish  et  al  (1971)  placed  several  deep  sea  bottom  pressure 
gauges  around  the  predicted  amphidrome.  The  record  from  one  gauge  (JOSIE  II 
■it  J.°N,  li5tW)  seaward  of  the  predicted  amphidrome  showed  that  indeed  there 
was  a  phase  difference  of  160°  between  this  location  and  measurements  on  the 
opposite  side  of  the  amphidrome.  However,  the  tide  constants  for  JOSlf.  II 
differ  from  that  predicted  by  Munk  cl  al  for  this  location  by  14  cm  in  amplitude 
and  18°  in  phase. 

The  solution  by  Pekeris  and  Accad  (1970)  is  a  numerical  integration  of  the 
Laplace  tidal  equat  ions  for  the  globe,  with  no  earth  tides,  but  with  tidal 
dissipation  parameterized  as  linearized  bottom  stress  in  shallow  water.  They 
i  i'd  a  fairly  detailed  model  of  the  ocean  bottom  topography,  but  did  not  use 
:  <!  'bsi- rvat  ions  of  tides  as  boundary  conditions.  They  found  that  the 
■resisted  tide  is  very  sensitive  to  changes  in  the  coastline  conf igurat ion . 


This  is  evident  in  Figure  2,  from  the  difference  in  amplitude  (14  cm)  and  phase 
(107°)  between  solutions  8  and  9,  where  the  coast  lines  are  modeled  on  a  1°  and 
2°  grid  respectively. 

Hendershott  [1972]  modified  Laplace's  tidal  equations  to  account  for  ocean 
bottom  loading  and  sel f-at t ract ion  as  well  as  the  astronomical  forces  in  the 
yielding  of  the  ocean  bottom  due  to  earth  tides.  Dissipation  was  parameterized 
as  a  near  coastal  process. 

Luther  and  Wunsch  (1975)  constructed  updated  tidal  charts  for  the  Pacific 
using  a  senl-empir leal  technique  similar  in  concept  to  those  of  Harris  and 
Dietrich  but  with  more  recent  tidal  observations. 

Estes  [197  7]  extended  the  approach  of  7-ahel  [1970]  to  the  0^ ,  ^  ,  K, ,  S.,  , 

and  N,  components  as  well  as  the  M.,  component.  Self-attraction  and  ocean  bottom 
1 oad in ,  were  also  Included  following  Hendershott  [1972].  These  effects  on  the 
ocean  tide  determination  can  be  quantified  by  comparing  determinations  12  and  13 
In  Figure  2.  A  shift  in  amplitude  and  phase  of  6  cm  and  27°  respectively  is 
observed  die  to  sc  1 f-at t ract ion  and  ocean  bottom  loading.  This  is  similar  in 
i  .  itude  t  the  change  in  the  determinations  14  and  15  by  Cordeev  et  al  [  1977  ] 
(3  cm  and  lfc°).  These  latter  investigators  also  used  self-attraction  and  ocean 
bott.v.  loading  following  Hendorshot t ,  but  the  treatment  of  dissipation  was 
similar  to  that  of  Pckcris  and  Accad,  i.c.,  bottom  stress  linear  with  velocity. 

The  final  technique  to  be  considered  is  that  of  Kuo  and  Jachens  [1977a] 
vh..  used  t'o  variation  of  gravity  on  land  to  deduce  the  amplitude  and  phase  of 
the  ocean  tide.  It  has  long  been  known  that  ocean  bottom  loading  affects 
surface  gravity  measurement s ,  even  at  intercontinental  distances.  By  a  linear 
programming  inversion  technique,  components  of  the  tides  can  be  mappel!  from 
t ldal  gravity  signals  at  a  network  of  terrestrial  stations.  The  results  of 
t  1-.  technique  ar<-  somewhat  sensitive  to  the  distribution  of  these  stations 
!  •  t hi  tide  r.»det  assumed  for  oceans  outside  the  region  to  be  modeled, 

:  u  . ic Id  good  agreement  with  independent  bottom  pressure  gauge 

t  is  .  renont  s. 


-9- 


With  the  exception  of  the  clot  crmin.it  Ion  number  lb,  which  i;,  1  rom  a  bott. 
pressure  gauge  located  at  the  compar Ison  point,  all  other  determlnat Ions  shown 
in  Figure  2  were  made  from  observations  which  are  considerably  distant  from  the 
comparison  point.  The  wide  spread  in  amplitudes  and  phases  for  this  point 
arc  due  primarily  to  differences  in  techniques  used  in  these  determinations. 

It  is  interesting  to  note  that  the  earliest  and  the  latest  determinations 
(Harris,  [1904],  and  Kuo  and  Jachens  [  197  7  )  are  only  14  degrees  different 
In  p!u»:;e.  One  might  conclude  that  tide  determinations  themselves  have  a  peri¬ 
od!  itv  of  70  years!  There  seems  to  he  a  slight  preference  among  the  determin¬ 
ations  for  the  27o'  phase  angle,  in  agreement  with  the  .IOSIF  II  data  hut  it  is 
hard  to  rule  out  determlnat Ions  2  through  5,  9  1  10,  and  14  &  15  as  being  invalid. 
:i  has  been  speculated  (Hunk,  et  al  (1970))  that  the  use  o!  coastal  tide 
observations  for  deep  ocean  tide  determination  Is  questionable-  because  of 
the  severe  local  anomalies  that  Can  occur  In  coastal  environments.  It  Is  a 
fact  that  most  coastal  tide  gauges  are  located  in  estuaries  or  bays  that  are 
protected  from  the  waves  and  winds  of  the  open  ocean.  They  are  also  protected 
from  the  open  ocean  tide  by  up  to  60  degrei s  of  phase  difference.  A  similar 
situation  exists  for  tide  gauges  on  islands.  They  arc  usually  placed  Inside 
lagoons  or  in  narrow  straits  between  islands  for  local  navigation  purposes. 

These  observations  can  have  30  to  60  degrees  phase  difference  from  the  tide 
in  ih.  adjacent  open  ocean  (Gallagher  et  al  (19711).  But  those  tide  determin¬ 
ation:.  reviewed  above  which  depend  on  these  observations  as  boundary  conditions 
.!*>w  no  distinctive  distribution  of  amplitude  or  phase  which  night  distinguish 
the  as  a  cla.s  1 roa  the  other  types  of  solutions.  Neither  do  t h<  empirical 
solution-.  .  ppear  better  or  worse  than  the  numerical  solutions.  There  appeal 
t  be  no  preferred  technique  among  the  tide-  determinations  examined  here. 

7  ie i r  variance  appears  to  arise  from  the  near  resonance  ol  the  tide  In  the 
-rthi  ist  Pacific  responding  to  small  differences  in  treatment  of  dissipation, 
and  i  .  probably  most  sensitive  to  the  paranct e r lz at  Ion  of  coast  lines. 

■  ;.li  i  -hott.  [1*7))  reco-mcnds  that  small  dot  ills  such  as  accurate  models 

c  .it  eut.tl  ielf  and  iland  scattering  and  el  if  ( rac  t  ie>n  he  developed 
re  proe.»  :tn,;  to  develop  deep  water  tide  riodcls  using  the  techniques 


.1  .•  r  i  bed  above . 


The  situation  la  ouch  the  same  ior  the  ten  1 -diurnal  tide  components 
(Figure  3).  The  number..  in  this  and  subsequent  figures  refer  to  the  legend 
In  Figure  2.  Again  a  wide  spread  in  amplitude  and  phase  is  evident,  with  no 
preferred  values.  Perhaps  the  S,  component  Is  also  near  resonance  in  the 
Northeast  Pacific.  The  other  t  Ide  components  examined  ;eem  to  show  preferred 
values  for  amplitude  and  phase,  although  this  is  Inconclusive  because  of  the 
smaller  number  of  determinations.  The  N,  component  (Fiore  4),  is  determined 
by  only  2  investigators,  und  they  report  values  of  4  rr  *  1  <  m  and  2f>bl  ■*  ]' 
(or  ,  amplitude  and  pfn.se.  The  diurnal  tide  components,  K,  and  Oj  (Figures 
5  and  b  respectively)  show  a  marled  preference  for  ?40‘  in  phase  and  29  cm 
and  1 b  cn,  respectively,  for  amplitude.  Presumably ,  the  Northeast  Pacliic 
Ocean  basin  is  not  resonant  for  these  components  of  the  tide. 

It  should  b»  noted  that  most  ol  the  e  tide-  parameters  have  been  extracted 
!to  publlsht  tide  charts  of  co-tidal  and  co-range  lines.  These  charts  have 
vary  In,  resolution  capability,  and  interpolated  values  of  amplitude  and  ptiase 
..te  !  .-ar  in,  accuracy.  A  conservative  estimate  of  tie  accuracy  of  these 
pal  u-.elel:.  would  be  *.  N  cm  (or  amplitude  an  I  t  S  for  phase  angle.  In  other 
word  ,  t  hi  accuracy  of  reading  these  chart.-,  usually  exceeds  the  tccuracy  of 
t  h>  charts  themselves,  especially  In  the  ca  .  of  the  t  ide. 

(liven  this  state-of-the-art  in  tide  determinations  in  the  Northeast 
Tacit  lc,  altimeter  tide  determinations  can  only  improve  our  knowledge  of  the 
X,  and  S,  components.  The  best  tests  for  evaluating  the-  altimeter  technique 
appear  to  be  comparison  with  K,,  Oj  and  perhaps  N,  determination-  iiscussco 
above . 
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ALTIMETER  DATA  PRO* ESSING 


The  Northeast  Pacific  region  is  traversed  by  about  200  CEOS-3  altimeter 
data  passes  (see  Figure  1),  and  about  400  Seasat  altimeter  data  passes.  Of 
these  data,  184  CEOS-3  passeB  have  already  been  processed,  and  the  data, 
reduced  to  sea  heights  and  corrected  for  orbit  errors,  is  available  from 
Goddard  Space  Flight  Center. 

To  appreciate  the  key  role  of  orbit  error  correction  In  the  recovery 
of  tidal  Information,  It  Is  helpful  to  review  the  altimeter  measurement 
geometry  shown  In  Figure  7.  The  altimeter  measures  the  distance  h  between 
the  Instantaneous  sea  surface  and  the  satellite  along  a  line  normal  to  the 
mean  sea  surface  at  a  given  time  t.  If  the  geocentric  distance  of  the  satel¬ 
lite,  r  ,  is  known  accurately  at  time  t,  then  the  geocentric  distance  to  the 
s 

Instantaneous  sea  surface  can  be  calculated.  In  particular,  the  sea  surface 
height,  S,  above  the  reference  ellipsoid  can  be  calculated. 

However,  if  there  Is  some  orbit  error,  l.e.,  rg  Is  not  known  accurately, 
then  the  corresponding  sea  surface  height  S  will  be  in  error  and  an  erroneous 
ocean  tide  height  will  result.  While  It  Is  theoretically  possible  (Zetlcr 
and  Maul  (1971})  to  recover  the  tide  parameters  In  the  presence  of  orbit 
errors.  Maul  and  Yanaway  (1978)  found  that  In  practice,  orbit  errors  were 
either  too  large  or  not  random  enough  to  permit  successful  tide  parameter 
recovery.  Other  investigators  (Won  and  Miller  (1978)  have  also  cited  exces¬ 
sive  orbit  errors  as  a  contributing  factor  In  failure  to  achieve  good  tide 
parameter  recovery  from  altimetry. 

In  this  investigation,  a  unique  process  is  used  to  remove  orbit  terrors 
from  the  altimeter  data.  As  it  Is  not  possible  to  know  accurately  the  geo- 
c entile  satellite  distance  rg  for  each  satellite  pass,  this  quantity  Is 
replaced  by  a  fictitious  quantity  which  will  result  In  the  minimum  deviation 
!  the  calculated  sea  surface  height  S  from  the  modeled  mean  sea  surface 
eight  N,  over  the  length  of  the  pass.  The  mean  sea  surface  (or  geoid)  height, 
,  is  calculated  along  the  satellite  subtrack,  N(t),  from  a  geopotential 
oc.,-1  such  as  GEM-7  (Wagner,  et  al  (1976)).  While  N  and  S  are  not  generally 
co-1 Inear,  to  a  good  approximation  they  may  be  assumed  so.  The  error  In  this 
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assumption  Is  In  the  sub-centimeter  range.  The  fictitious  geocentric  satellite 
distance  r  may  be  parameterized  as  a  low  order  polynomial  in  t,  c.g. , 

r#(t)  a  a  +  bt  +  ct^ 

over  the  relatively  short  length  of  the  satellite  pass.  This  is  a  good  approxi¬ 
mation  since  only  snail  fractions  of  a  complete  satellite  orbit  are  to  be  fitted 
In  each  pass.  The  polynomial  parameters,  a,  b,  c,  are  then  adjusted  In  a  least 
squares  process  to  minimize  the  difference  between  the  calculated  sea  surface 
height,  S(t) 


S ( t )  •  r  (t)  -  h(t)  -  r  (t) 

s  e 

and  the  calculated  geold  height,  N(t),  along  the  length  of  the  pass.  Since 
each  pass  Is  thus  fitted  independently  to  the  some  model  mean  sea  level  surface, 
all  orbit  errors  and  any  other  phenomena  which  manifests  Itself  like  an  orbit 
error  (l.c.,  a  bias  or  tilt)  are  removed  and  the  sea  height  difference  repre¬ 
sents  the  true  difference  In  sea  surface  height  at  the  times  of  the  intersecting 
passes.  Each  pass  uses  an  independent  set  of  fictitious  polynomial  parameters, 
an.!  these  have  no  real  meaning  or  physical  function  whatever,  except  to  accom¬ 
plish  the  fitting  of  the  sea  surface  heights  to  the  model  mean  sea  level.  In 
this  process,  care  must  be  taken  In  choosing  the  length  of  the  pass  which  Is 
fitted.  In  a  very  short  pass,  even  the  ocean  tide  signal  appears  as  a  bias 
and  would  be  removed  In  this  process.  Therefore,  all  passes  arc  chosen  to 
be  longer  than  2500  km  (6  minutes),  which  Is  sufficient  to  preserve  the  tidal 
information  of  the  Northeast  Pacific  ocean  basin. 

In  actual  practice,  the  adjustment  and  fitting  of  the  sea  surface  height 
■ita  is  accomplished  in  a  general  purpose  orbit  determination  program,  using 
r«  soph  1 »t icat ed  calculation  of  r  ,  but  the  principal  is  the  same.  It  is 
.  tir-ited  that  this  process  results  In  sea  surface  heights  which  are  free  of 
r  it  error  to  wlthir  20  centimeters  (Brown  et  al  (1977)). 
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CROSSOVER  I"T\  PROCESSING 


Once  the  altimeter  data  are  converted  to  corrected  sea  surface  heights, 
it  Is  possible  to  parameterize  difference.?  in  sea  heights  at  crossovers  in 
terms  of  phase  differences  of  various  components  of  the  tides,  the  five 
largest  tide  components,  K^,  S^ ,  Oj ,  and  N^,  are  estimated  in  a  harmonic 

analysis  of  the  crossover  discrepancies.  Referring  to  Figure  6,  each  sea 
height  is  expressed  as  the  sum  of  contributions  of  all  5  components. 

5  5 

5(t)  -  l a ^  sin  w^t  +  cos  w^t]  ■  cosfw^t  -  ♦  )  (1) 

1-1  i-1 


.  rly,  to  solve  for  the  ten  unknown  parameters,  we  need  at  least  10  indepen- 
t  measurements  of  AS  at  a  particular  location.  This  is  difficult  to  achieve, 
. '  \  t  n  that  the  CF.OS-3  satellite  orbit  never  exactly  repeats.  However,  it  does 
cone  close  to  repeating  (within  1/2  degree)  every  fortnight.  Thus,  it  is  not 
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difficult  to  collect  a  sufficient  number  of  crossovers  within  a  few  decrees 
of  the  desired  locution  (see  Figure  9)-  But  one  must  assume  th  t  the  tide 
amplitude  and  phase  are  constant  over  the  area  spanned  by  this  collection 
of  crossovers. 

Further  assumptions  are  necessitated  by  the  fact  that  one  does  not  collect 
a  sufficient  number  of  crossovers  Instantaneously.  Due  to  the  power  limitations 
of  the  CEOS-3  satellite,  the  altimeter  did  not  operate  continuously.  Therefore, 
there  nay  be  a  time  lapse  of  a  year  or  more  between  1 7  and  t  j  at  a  given  cross¬ 
over  (see  Figure  9).  Thus  one  must  assume  that  there  are  no  secular  effects 
on  the  tides  which  are  significant  in  the  lime  span  needed  to  collect  sufficient 
crossovers. 

Such  long  phase  delays  make  accurate  knowledge  of  the  tidal  component 
periods  more  important  also.  An  error  of  about  one  part  in  10*  in  the  period 
of  a  semi-diurnal  component  is  sufficient  to  cause  a  2.5°  error  in  ph-.se  angle 
after  a  year.  Tide  component  periods  used  In  this  Investigation  nre  accurate 

g 

to  better  tlun  1  part  In  10  . 

Other  factors  that  change  with  time  must  also  be  carefully  mon ti< .red . 

For  example,  while  the  altimeter  is  largely  insensitive  to  sea  state  (wind 
waves,  swells,  etc.),  a  significant  wave  height  (SVH)  of  10  m  or  more  can 
Cause  J.5  a  to  l  a  error  in  altitude  readings  (Hofmelster,  et  al  (1976]). 

It  is  therefore  important  that  the  sea  state  conditions  hi  rather  calm 
(SWH  <  3m)  for  good  tide  determinations. 

The  OFOS-3  and  SEASAT-1  altimeter  data  d  i  st  r  iluit  ions  have  been  examined 
t  Vt  ermine  the  number  and  location  of  potential  tide  solution  locations 
(n  the  Northeast  Pacific.  Locations  at  which  concentrations  of  twenty  or 
,.,oi  crc:,„.;vcrs  exist  within  a  small  local  area  (less  than  about  5  by  5° 

«  i)  from  passes  longer  than  2500  km  have  been  Identified  for  both  satellites. 

.  ••(>•-!  lx  tide  solution  areas  which  have  been  identified  from  98  0F.0S-3 
a  .  art  shown  bv  the  circles  In  Figure  10.  The  solid  triangles  in  Figure  lr 
.  .sent  the  solution  areas  of  highest  crossover  density  from  the  SFASAT-1 
,ilt.  eter  data.  Since  the  SEASAT  satellite  was  placed  in  a  3-day  repeating 
r! it  nv.it  the  end  of  its  lifetime,  crossover  concentrations  in  excess  of 
I0f  per  6  min.  by  6  min.  area  were  generated  at  these  locations.  Many  other 

solution  locat  Ions  which  have  fewer  crossovers  could  bo  Identified  from  the 
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S1ASAT-1  data  distribution,  but  these  37  will  provide  potentially  tie  highest 
quality  solutions.  The  total  of  63  tide  solution  locations  identified  in 
Figure  10  provide  for  a  tide  deterrainat ion  roughly  every  5  degrees. 

The  determination  of  tide  parameters,  and  b^,  at  a  given  location 
is  a  staple  least  squares  adjustment  to  minimize  the  crossover  discrepancies 
ASj .  Normal  equations  are  formed  for  each  of  the  solution  locations.  If  we 
rewrite  equation  2  in  terns  of  matrices,  we  have 

AS  •  (A)  X 


where 


and 


IJ 


2  1 
is  sin  j t  j  —  Rin  « 1 1.  j , 


2  1 

or,  cos  wjlt  "  cos  wjlj 


The  resulting  normal  equation 

|a]t  AS  -  (a]t  (A)  X 


ha  .  as  it  s  solut ion 


PRELIMINARY  TILE  SOLUTIONS 


Tide  determinations  have  been  made  to  date  at  eight  of  the  26  t.EOS-3 
solution  locations  saown  in  Figure  10.  These  solutions  are  of  a  preliminary 
nature,  but  serve  to  point  out  some  problem  areas  as  well  as  some  encouraging 
new  approaches  to  be  investigated  in  a  subsequent  study. 

Table  1  shows  the  results  of  tide  solutions  at  several  of  these  locations 
for  up  to  5  tide  components.  All  of  the  solutions  are  standard  unconstrained, 
unweighted  least  squares  solutions.  The  first  two  solutions,  TMA  and  TMB , 
illustrate  a  phenomenon  which  nay  be  called  dilution  of  the  strength  of  the 
solution.  Both  solutions  result  in  a  reduction  in  the  RMS  of  crossover 
residuals,  but  at  the  expense  of  ridiculously  large  tide  amplitudes.  The 
determinant  values  are  also  relatively  small  (even  negative  in  the  case  of 
T>3)  which  further  Indicates  a  weak  or  Ill-conditioned  solution.  This  situa¬ 
tion  cones  about  when  the  number  of  unknowns  (10)  approaches  the  number  of 
obse rvat  Ions  (27  and  25).  It  is  Improved  as  the  number  of  unknowns  is  decreased. 
For  example,  note  the  decrease  In  amplitude  and  Increase  In  determinant  value 
as  the  number  of  conponcnis  solved  for  In  location  2  Is  reduced  (solutions  TM 2, 
:M5,  TM52 ,  and  TMB2).  The  value  for  M.,  achieved  in  TMB2  is  in  fair  agreement 
w'th  bottom  pressure  gauge  results  reported  by  Munk  et  nl  (1970j  for  the  nearby 
Fllloux  gauge  (.19  n  amplitude  and  107.14°  Creenwich  phase  angle). 

This  approach  of  reducing  the  unknows  is  not  generally  recommended  because 
■a*  everal  frequencies  represented  In  the  5  tide  components  are  not  distinctive 
>.  that  the  parameters  of  each  component  are  uncorrclat ed .  Thus,  If  a  given 
,  nent,  say  N ,  is  not  solved  for.  Its  contribution  to  the  crossover  discrcpan- 
..  1.  absorbed  erroneously  by  another  component  whose  parameters  have  a  high 

.!  -c  rrelation  with  N,.  This  is  seen  in  the  dramatic  change  of  parameter 
val  ..  caused  by  deletion  of  the  0^  component  in  solutions  TT12  and  TM5.  In 
■4-1  t  :  'ti  X)'.2 ,  the  correlation  between  the  0j  parameters  and  other  parameters 
ranges  from  .81  to  .91. 

The  addition  of  further  observations  also  strengthens  the  solution,  as 
'vi  in  mo  1  ii t  Ion  ADD,  but  not  as  dramatically  as  reducing  the  unknowns.  The 
r-.il  equit Ions  of  solutions  TMA  and  TMB  were  added  and  solved  to  produce  ATD. 

!  I Is  equivalent  to  a  single  solution  of  5?  cross-over  discrepancies.  Note 
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PRE!  1  MI  NARY  TIDE  SCLUTIOriS 


that  the  determinant  value  increased  and  the  tide  component  amplitudes  decreased 
to  more  reasonable  levels.  This  is  a  preferred  technique  for  improving  the 
qu~’lty  of  the  tide  solutions.  In  this  case,  the  addition  of  observat ions  is 
mi; leading,  since  they  are  not  from  the  name  immediate  area.  In  ADD,  the  cross- 
ovi  locations  range  from  22.5°  to  30.5°N  and  224°E  to  ?29°E.  When  SEASAT  data 
an  processed,  the  crossovers  available  for  a  given  location  are  expected  to 
increase  by  an  order  of  magnitude,  with  a  corresponding  increase  In  quality  of 
the  solution. 

Another  preferred  technique  for  improving  solutions  is  to  imj rove  the 
quality  of  the  observations.  In  a  few  solutions,  it  was  noticed  that  the 
crossover  discrepancy  values  for  a  given  pass  were  systematically  different 
from  those  of  other  passes,  e.g.,  all  crossovers  associated  with  a  given  pass 
night  have  a  significant  non-rero  mean.  This  could  happen  if  orbit  error 
removal  was  incoaplcte  for  some  reason.  When  these  crossovers  are  removed 
from  the  solution,  the  solution  values  became  more  reasonable,  despite  the 
reduced  number  of  observations.  For  example,  solution  TM3  for  location  3 
had  bO  crossovers,  hut  those  of  3  passes  appeared  to  be  syst cmat it al ly  offset 
: ror:  the  others.  When  these  crossovers  were  removed,  as  in  solution  TM3L, 
the  amplitude  of  the  term  dropped  from  2.1  to  .3  meters.  The  1-MS  of  cross¬ 
over  discrepancies  also  improved  more  in  TK3L  than  in  TM3,  despite  the  loss 
of  30  observations.  Apparently,  biases  in  the  crossover  discrepancies  alias 
significantly  with  tide  parameters. 

Ideally,  one  would  like  to  be  able  to  correct  the  biases  in  faulty  passes 
while  still  retaining  the  tidal  information.  This  may  be  e'foctet!  by  adding  pass 
la.  terms  to  the  observation  equation.  In  this  case,  equation  (?)  becomes 

5 

•  0,  -  -*•  £  |"at(8in  w^t.,  “  sin  w^)  +  kj(cos  witi  “  cos  wiCi^J  ^') 

'.'or  al  matrices  formed  from  these  equations  would  have,  in  addition  to  the 
10  tide  parameters,  one  pass  bias  for  each  pass  (one  pass  bias  rust  be  held 
fixed  to  avoid  a  singular  matrix).  This  has  the  effect  of  diluting  the 
strength  of  the  solution  if  all  parameters  are  solved  simultaneously,  hut 
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It  la  feasible  to  solve  for  the  biases  Independently  and  correct  the  cr >ssover 

discrepancies  before  executing  the  tidal  solution.  This  was  attempted  In 

solution  T04B.  Comparing  with  solution  T04 ,  before  pass  bias  correction,  we 

see  that  Che  most  significant  effect  Is  the  reduction  of  the  amplitude  of  the 

R,  component.  Not  much  phase  angle  change  occurs  at  all. 

«• 

Finally,  the  solution  for  location  8  Is  presented  to  indicate  that  there 

are  data  which  yield  well-behaved  solutions  without  any  remedial  tteatment. 

This  Is  a  case  where  the  input  data  are  not  biased  (Initial  RMS  of  crossover 

is  small,  +  1.52  a),  and  there  are  sufficient  crossovers  (34)  to  give  a  strong 

12 

solution  (determinant  is  rather  large  at  1.1  x  10  ). 

It  Is  too  soon  to  begin  comparison  of  the  recovered  tide  parameter  values 
with  those  of  other  Investigators,  but  the  preliminary  results  are  encouraging. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Knowledge  of  tides  In  Che  Northeast  Pacific  Ocean  is  fairly  uncertain  at 
present,  but  ihe  advent  of  satellite  altimetry  data  from  the  CEOS-3  and  SEASAT 
satellites  promise  the  possibility  of  direct  measurements  of  tides  in  this  area. 
Kror  these  data  one  can  determine  tide  Constants  on  a  5°  by  5°  grid  over  all  the 
deep  ocean.  Analysis  of  preliminary  altimeter  tide  solutions  have  underscored 
the  necessity  for  complete  removal  of  all  orbit  type  bias  errors  from  the  data 
before  attempting  to  solve  for  tides.  Fortunately,  techniques  have  been  developed 
for  editing  or  correcting  biased  data  and  obtaining  good  quality  tidal  solutions. 
The  future  processing  of  SEASAT  data  will  yield  much  higher  quality  solutions 
owing  to  the  high  density  of  crossovers  which  accrues  from  the  SEASAT  repeating 
orbit,  and  the  more  accurate  SEASAT  altimeter  data. 

In  the  near  future,  the  remaining  CEOS-3  crossover  data  will  be  processed 
and  corrected  for  biases  and  solutions  will  be  achieved  for  all  26  solution 
areas.  Processing  of  SEASAT  altimeter  data  is  expected  to  begin  in  November 
of  1979.  Comparison  and  validation  of  various  tide  determinations  will  follow 
the  construction  of  tide  solutions  from  the  SEASAT  data. 
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*  Tides  in  the  deep  ocean  can  be  dctirmincd  directly  from  satellite 
altimetry,  completely  Independent  of  assumptions  about  earth  tides, 
bottom  topography  and  coastal  geometry  and  thus  free  of  the  uncertainties 
which  plague  numerical  tide  models.  Existing  tide  models  differ  by  1  meter 
or  norc  in  the  value  of  sea  surface  height  in  the  deep  ocean  at  a  given 
, lace  and  time.  This  uncertainty  is  a  formidable  bstacle  to  determination 
.1  a  precise  marine  gcolds  t rom  satellite  altimetry.  By  harmonic  analysis 
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of  the  temporal  ciianges  In  nltimeter  measurements  at  satellite 
subtrack  crossover  points.  It  Is  possible  to  so’ve  for  the 
amplitude  and  phase  of  harmonic  tidal  components.  However, 
care  must  be  exercised  In  the  removal  of  satellite  orbit  errors, 
and  In  the  selection  of  crossovers  for  sufficient  observability  of 
the  phase  angle  of  the  harmonic  tidal  component.  Preliminary  tidal 
solutions  In  the  Gulf  of  Alaska  using  the  relatively  spcrac  CEOS-3 
altimeter  data  distribution  show  generally  good  agreement  (20  cm  in 
amplitude  and  25  degrees  in  phase)  with  deep  ocean  bottom  pressure 
gauge  measurements  and  stabllsh  the  feasibility  of  this  technique. 
SEASAT  altimeter  data  yields  a  much  greater  density  of  crossovers 
(400  plus  per  1/2*  by  1/2P  area),  making  possible  much  better 
separation  of  Individual  harmonic  components  within  the  semi¬ 
diurnal  and  diurnal  families.  .  .Five  degree  resolution  on  tidal 
charts  of  co-range  and  co-phase  lines  in  the  Gulf  of  Alaska  are 
possible  with  the  SEASAT  data  distribution. 
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